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Preparation and characterization of arbitrary states of four-dimensional qudits based 

on biphotons 
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We report interferometric schemes to prepare arbitrary states of four-dimensional qudits 
(ququarts) based on biphoton states of ultrafast-pumped frequency-nondegenerate spontaneous 
parametric down-conversion. Preparation and tomographic characterization of a few examples of 
general single-ququart states, a pure state, a completely mixed state, and a partially-mixed state, 
are experimentally demonstrated. 

PACS numbers: 03.67.-a,03.65.Wj, 42.50.Dv, 42.50.-p 



INTRODUCTION 

In quantum information, a two-dimensional quantum 
system is often used as a carrier of the basic information 
unit, the two-level quantum state or the quantum bit 
(qubit). Most quantum computing and quantum com- 
munication protocols are based on preparation, manipu- 
lation, entanglement, distribution, and measurement of 
multiple qubits [Ij]. 

Recently, D-level quantum states [D > 2) or qudits 
have attracted a lot of attention in the context of quan- 
tum communication and experimental tests of quantum 
mechanics 0, S i, i, S 0, H, Hi- Experimentally, the 
physical carrier of the qudit can be any D-dimensional 
quantum systems. Meaningful applications of qudits in 
quantum information, however, will only be possible if 
the qudit is encoded in a D-dimensional physical degree 
of freedom which is easy to handle experimentally. Fur- 
thermore, it should be possible to entangle multiple in- 
dividual qudits in a scalable manner. 

In photonic quantum information research, an internal 
or an external degree of freedom of a photon is used to 
encode the intended quantum state. For a qubit, the po- 
larization state of a photon is often the obvious choice 
although it is possible to choose other degrees of free- 
dom. To encode a qudit, it is necessary to choose a 
multi-dimesional degree of freedom of a single-photon, 
such as, the angular momentum, transverse momentum- 
position, time of arrival, etc 10|, luJ, [l^, llJ, llJ, uM llSf . 
These single-photon multi-dimensional photonic degrees 
of freedom, however, are experimentally difficult to ma- 
nipulate and there are no scalable schemes to generate 
multi-qudit entangled states. 

The single qudit, however, does not need to be en- 
coded in a single-particle quantum state. In fact, prepa- 
ration and tomographic characterization of a pure state 
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qutrit (three-dimensional quantum state), see Refs. 17|, 
ISl . Il9|, and a pure state ququart (four-dimensional 
quantum state), see Refs. [20|, l2l|, have been demon- 
strated recently using the biphoton polarization states 
of frequency-degenerate and frequency-nondegenerate 
spontaneous parametric down-conversion, respectively. 
In other words, a pair of photons can be used as a carrier 
of three- or four-dimensional quantum states. 

Especially, the ququart based on the biphoton po- 
larization state of frequency-nondegenerate spontaneous 
parametric down-conversion (SPDC) exhibits a few prop- 
erties which are important for applications in quantum 
information research: First, all the ququart basis states 
can be accessed using only linear optical elements (phase 
plates) 20,l21[. Second, it is possible to prepare a multi- 
ququart entangled state starting from multiple individual 
biphoton ququarts, linear optical elements (beam split- 
ters), and post-selection measurement [22|, |23|. So far, 
only pure state biphoton ququarts have been demon- 
strated. It is, therefore, of interest and importance 
to learn how to prepare a general state of a biphoton 
ququart [24 1. 

In this paper, we report experimental studies on 
preparation of general states of biphoton ququarts using 
ultrafast-pumped frequency-nondegenerate spontaneous 
parametric down-conversion. Methods for preparation 
and tomographic characterization of some examples of 
arbitrary ququart states, i.e., pure, mixed, and partially- 
mixed states of a biphoton ququart, are experimentally 
demonstrated. We also discuss a couple of alternative 
experimental schemes which allow to generate arbitrary 
biphoton-based ququart states. 



BIPHOTON QUQUART 

Let us start with a brief introduction to the biphoton 
ququart. In coUinear frequency-nondegenrate SPDC, a 
higher energy pump photon is occasionally split into a 
pair of co-propagating lower energy photons (signal- idler) 



of different frequencies, 

1 1 



A. 



Ai A2 



(1) 



where Ap is the pump wavelength and Ai (A2) is the wave- 
length of the signal (idler) photon. Since each photon of 
the pair can be horizontally or vertically polarized, the 
following set of biphoton polarization basis states can be 
defined [23, [m. Si [13, 

\Hx,,Hx,) = \0), \Hx,,Vx,) = \l), 

The polarization state of the photon pair born in the pro- 
cess of collinear frequency-nondegenerate SPDC, there- 
fore, represents a four-dimensional quantum state or a 
ququart. Since the states shown in eq. ([2]) are orthonor- 
mal to each other and form a complete basis for a four- 
dimensional Hilbert space, they form the computational 
basis for the biphoton ququart. 

Experimentally, the computational basis states for the 
biphoton ququart can be generated using type-I SPDC 
(|0) and |3)) and type-II SPDC (|1) and |2)). To prepare 
an arbitrary superposition state of a ququart. 



IV) =co|0)+ci|l)+C2|2)+C3|3), 



(2) 



where c; = jcjle**^' are the complex probability ampli- 
tudes which satisfy J2i=o l*^'!^ ~ ^' '^^'^ would need to 
coherently combine two type-I SPDC and two type-II 
SPDC sources and to control four complex amplitudes 
Ci,C2,C3, and C4. The use of four SPDC sources, how- 
ever, turns out to be unnecessary since the degree of 
polarization of the biphoton state of collinear frequency- 
nondegenerate SPDC is not invariant under SU(2) trans- 
formations [25|]. It is, therefore, possible to use linear 
optical elements (wave plates) to transform one bipho- 
ton ququart basis state into any other one as recently 
demonstrated in Ref. [20|, |2]J. As a result, it is possi- 
ble to prepare an arbitrary superposition state (i.e., pure 
state) shown in eq. ^ using less than four SPDC sources 
and we discuss several such schemes in this paper. 

The superposition state shown in eq. ^, however, 
is not the most general quantum state for a ququart. 
To properly consider mixedness, which comes from con- 
trolled or uncontrolled quantum distinguishability among 
the basis states shown in eq. ([2]), the ququart state 
should be expressed as a 4 x 4 density matrix p to de- 
scribe general single-ququart states: completely mixed 
states, partially-mixed states, and the pure state shown 
in eq. p. 

It is interesting to note that, since the ququart state 
under consideration is in fact formed with two polariza- 
tion qubits, it is possible to apply the state classifica- 
tion method based on two-qubit entanglement of forma- 
tion. This method is based on the quantity C, concur- 
rence, which is a measure of two-qubit entanglement [26l | . 
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FIG. 1: Schematic of the experiment. Ququart preparation: 
The collinear frequency-nondegenerate biphoton polarization 
state is transformed by the wave plate (WP) to prepare an ar- 
bitrary four- dimensional quantum state, a biphoton ququart. 
The theoretically expected quantum state of the ququart, 
Ptheory, Can be calculated by using the information on the 
experimental settings. Ququart tomography: The biphoton 
ququart is tomographically characterized to obtain the ex- 
perimentally reconstructed density matrix pexp- See text for 
details. 



For the single biphoton ququart (i.e., biphoton two-qubit 
state) in eq. ([2]), it is easy to show that there is no two- 
qubit entanglement (C = 0, i.e., the two-qubit state is 
separable) if the amplitudes satisfy the relation C0C3 — 
C1C2. Otherwise, the state is an entangled two-qubit state 
with concurrence C = 2|coC3 — C1C2I > 0. (The two-qubit 
state is non-maximally entangled if < C < 1.) 



QUQUART PREPARATION 

In order to prepare a general single-ququart state, in- 
cluding mixed, partially-mixed, and pure states, it is nec- 
essary to introduce quantum distinguishability among 
the biphoton ququart basis states defined in eq. (U) in 
a controllable manner. Moreover, to establish a confi- 
dence bound and to find a reliable method of preparing 
an intended ququart state, it is required to compare the 
experimentally reconstructed and the theoretically ex- 
pected ququart density matrices. Our experiment, there- 
fore, deals with both the ququart preparation as well as 
the experimental reconstruction of the ququart density 
matrices using the quantum state tomography. 

The schematic of the experimental setup is shown in 
Fig. [TJ In the single- crystal scheme, a type-I BBO crys- 
tal, with its optic axis horizontally oriented, is pumped 
by a train of ultrafast pulses and this scheme is used for 
preparing pure state ququarts. For preparing mixed and 
partially-mixed states, we used the double-crystal scheme 
in which two orthogo nally oriented type-I BBO crystals 
are placed in tandem [27|. The pump laser had roughly 
100 fs pulse width and centered at 390 nm. The polar- 
ization of the pump was controlled by using a half-wave 
plate (Pump HWP). The BBO crystals used on this ex- 
periment are 3 mm thick. The SPDC signal-idler photon 
pair generated at the crystals propagates coUinearly with 
the pump and has different frequencies: the signal pho- 



ton is centered at 823.5 nm (Ai) and the idler photon is 
centered at 740.8 nm (A2). 

In the single-crystal scheme, the biphoton polarization 
state generated at the crystal is IVaiiVa^) which corre- 
sponds to the ququart basis state |3) in eq. ^. This 
basis state is then unitarily transformed into a superpo- 
sition state shown in eq. ^ to prepare an arbitrary pure 
state by using a zero-order wave plate (WP). This is due 
to the fact that the signal and the idler photons, even 
though they have the same polarization, acquire differ- 
ent phase shift due to the large difference in the wave- 
length. As a result, the polarization states of the signal 
and the idler photons evolve differently, but predictably, 
for a specific angle of WP: \Vx,) -* a\Hx,) + l3\Vx,) 
and \Vx,) -^ "f\Hx^) + 6\Vx,) with \a\^ + \P\^ = 1 and 
ItP + I'^P = 1- The final result is the transformation of 
the initial ququart basis state into a superposition state 
in cq. ^ which is a pure state. The theoretically ex- 
pected density matrix ptheory for the prepared ququart 
state, which can be calculated accurately as the settings 
of WP is in our control, is then compared to the experi- 
mentally reconstructed density matrix pexp- 

The single-crystal scheme, therefore, allows us to pre- 
pare one of the four single-ququart basis states in eq. ^ 
initially. For the case just described, the single-ququart 
basis state |3) is in fact \V\^ , VX2), which is a factorizable 
two-qubit state with C — 0. Note now that the degree of 
entanglement (in this case, between the two polarization 
qubits) cannot be increased or decreased by local unitary 
transformations. The unitary transformation due to WP 
can turn the initial ququart state |3) in to a superposition 
state in the form of eq. ([2]). However, the resulting state 
can always be expressed as a separable two-qubit state 
of the form \Ax-^,B\2), where \A\^) and JBas) represent 
arbitrary polarization states of signal and idler photons, 
respectively. 

For preparation of mixed and partially- mixed ququart 
states, the double-crystal scheme described earlier is 
used. Since the two BBO crystals are orthogonally ori- 
ented, it is possible to excite two of the ququart ba- 
sis states shown in eq. ([2]): |0) (IHx^^Hx^)) and |3) 
{\Vx^,V\2)). The relative amplitudes between the two 
can be controlled by changing the pump polarization. 

To prepare a completely mixed state of |0) and |3), no 
further actions are required as the two amplitudes are 
already distinguishable in time due to the clock effect of 
the pump pulse [27.]. The resulting density matrix is, 
therefore, 



meo.,= (l-|)|0)(0| + ||3)(3|. 



(3) 



where the parameter x can be varied by changing the 
polarization of the pump. Adding two additional type- 
II BBO crystals, which are orthogonally oriented, will 
allow us to easily prepare a mixture of all four ququart 
basis states. This, however, is not necessary in principle 



since it is possible to transform a single basis state into 
a superposition of all basis states and then to introduce 
bircfringent/dichroic decoherence among the amplitudes. 
A more general state, between the completely mixed 
and the pure states, would exhibit some coherence among 
the four ququart basis states in eq. ^. In other words, 
the ququart density matrix has non-zero off-diagonal el- 
ements. Such states can be prepared by unitarily trans- 
forming the mixed state in eq. ([3]) using WP. As previ- 
ously discussed, a ququart basis state can be transformed 
into a superposition of all basis states linear optically us- 
ing WP [25[. By subjecting the mixed state in eq. ([3]) to 
unitary transformation using WP, it is possible to obtain 
the following partially-mixed state. 



Pth 



eory 



Pl\lpl){lpl\ +P2\^2){^2 



(4) 



where pi + P2 — ^ and |'0)i7 for example, is in the form 
of eq. (12). It is important to note that this unitary 
transformation process, however, does not actually de- 
crease entropy of the ququart state as we shall show in 
the next section. To decrease the entropy, it is neces- 
sary to erase the temporal distinguishability of the bipho- 
ton amplitudes born in the first and the second crystals, 
for example, by inserting a piece of thick compensating 
quartz crystal in the pump beam or in the path of the 
photon pair [27|. Introduction of controllable birefrin- 
gent/dichroic decoherence will transform eq. ^ into a 
more complex ququart state with increased entropy. 

From the two-qubit perspective, linear optical state 
transformation from eq. ([3]) to eq. (|4|) represents no in- 
crease in the degree of two-qubit entanglement as both 
states exhibit the two-qubit concurrence C — 0. This is 
closely related to the fact that the ququart state entropy 
remains the same for states in eq. ^ and eq. ^. To 
be more specific, the two-qubit concurrence C will in- 
crease if the single-ququart entropy is decreased by eras- 
ing the temporal distinguishing information present in 
eq. dS]) or in eq. ([4]). Complete erasure of the temporal 
distinguishability between the biphoton amplitudes from 
the first {\H\^,H\.2)) and the second (|Vai,Vx2)) type-I 
BBO crystal, see Fig. [1] will result a pure ququart state 
with the two-qubit concurrence C = 1. 

In this paper, we have experimentally demonstrated a 
pure, a mixed, and a partially mixed biphoton ququart 
states which are shown in eq. ^, eq. ([3]), and eq. ^, 
respectively. For the mixed and partially mixed states, 
the state entropy can be controlled by inserting a proper 
compensating crystal before or after the BBO crystals. 



QUQUART TOMOGRAPHY 

The prepared ququart state is characterized experi- 
mentally by performing quantum state tomography, a 
statistical method of reconstructing the quantum state 



density matrix based on a set of polarization projection 



measurement 

The experimental schematic for the ququart state to- 
mography is shown in Fig. [T] First, the photon pair 
that forms the biphoton ququart is split into two spatial 
modes by using a dichroic beam splitter (DBS), which 
transmits Ai = 823.5 nm and reflects A2 = 740.8 nm. 
Each photon, then, undergoes polarization state trans- 
formation with the use of a quarter wave plate (QWP) 
and a half wave plate (HWP). Finally, the polarization 
state projection is applied to each photon by using a po- 
larizer (POL) after which the photon is detected at the 
detector package, which consists of a spectral filter and a 
multi-mode fiber coupled single-photon counting module 
(SPCM). Since the ququart is made of a pair of pho- 
tons, the ququart detection is based on the coincidence 
counting rates of the two SPCM's that detects Ai and 
A2 photons with definite polarizations. The coincidence 
window used in this experiment was 5 ns. 

As noted in Ref. [20|, the ququart based on 
the collinear frequency-nondegenerate SPDC photon 
pair is mathematically equivalent to the non-coUinear 
frequency-degenerate SPDC photon pair which is often 
used in quantum information research. It is thus possi- 
ble to apply the two-qubit quantum tomography method 
described in Ref. [28| directly in this experiment to re- 
construct the single-ququart density matrix. Sixteen par- 
ticular joint biphoton polarization state measurements 
are, therefore, necessary to reconstruct a single-ququart 
density matrix for the identically prepared ensemble of 
biphoton ququarts. Table [J shows experimental settings 
of WP to perform the sixteen polarization projection 
measurements. The set of sixteen coincidence measure- 
ment outcomes n^ allow the linear tomographic recon- 
struction of the single-ququart density matrix. It is, 
however, possible that the mathematically reconstructed 
density matrix by the linear tomographic reconstruction 
might violate the physical properties of a density matrix. 
To avoid this problem, the maximum likelihood method 
was applied as follows [285. First, we generate a physical 
density matrix which satisfies normalization, Hermitic- 
ity, and positivity, as a function of sixteen variables. We 
then introduce the "likelihood function" which quanti- 
fies how good physical density matrix is in relation to 
the experimental data. Finally, using standard numeri- 
cal optimization techniques, we obtain the best estimate 
of the density matrix by maximizing the likelihood func- 
tion. The single-ququart density matrix initially recon- 
structed by the linear tomography is used as the seed for 
the iteration algorithm. 

The experimentally reconstructed ququart density ma- 
trix pexp is then compared to the theoretically expected 



density matrix ptheory which is calculated from the known 
values of the pump polarization and WP settings. The 
fidehty F = (Tr y/ ^/pTheoryPexp^/ Ptheory^ IS then calcu- 
lated to see how closely the two overlap and the state pu- 

TABLE I: QWP and HWP (fast axis) settings for ququart 
tomography. POL in Fig. [l] transmits vertical polarization. 

Here, \D) = {\H) + \V))/V2, \A) = {\H) - |V))/^/2, \D) = 
{\H} + \V))/V2, and \R) = {\H} + i\V})/V2. 
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rity is analyzed by calculating the state entropy defined 
as 5* = — X]fe=i ^k log4 Afe, where Aj, are the eigenvalues 
of the density matrix p. 



Pure state ququart 

To prepare a pure state ququart as in eq. ^ , we used a 
single type-I BBO crystal generating the initial ququart 
state |3), which in fact is a factorizable two-qubit state 
|^Ai,^A2)- This state is then transformed into a super- 
position state with the help of WP shown in Fig. [TJ In 
this experiment, a zero-order half wave plate designed at 
823.5 nm was used in place of WP. As a demonstration 
of pure state ququart preparation, we set the fast axis of 
the half- wave plate at 30° from the vertical axis. Since 
local unitary transformations do not change the degree of 
entanglement, the final states belong to the factorizable 
subset of two-qubit (ququart) states with C — Q. 

The theoretically expected ququart density matrix in 
this case is calculated to be, 



pure 
r theory ~~ 



( 0.5432 0.3136 + 0.1182i 0.3136 0.1811 + 0.0683i \ 

0.3136 -0.1182i 0.2068 0.1811 - 0.0683i 0.1194 

0.3136 0.1811 + 0.0683i 0.1811 0.1045 + 0.0394i 

V 0.1811 -0.0683i 0.1194 0.1045 - 0.0394i 0.0689 

I 



(5) 
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FIG. 2: Real part of the experimentally reconstructed (peip) 
and the theoretically expected (ptheory) density matrices for 
a pure state ququart. The calculated fidelity is F = 0.938. 



For this state, it is easy to see that Tr[p^] = 1. The 
entropy of the theoretical density matrix p^^g^^^ is calcu- 



lated to be 5 = 0, as it should for a pure state. 

To obtain the ququart density matrix, we have per- 
formed the sixteen projection measurement described in 
Table |T1 The coincidence data (for the accumulation time 
of 180 s) are ni = 6118, na = 1858, ng = 917, rn = 2943, 
715 = 1565, ne = 477, nr = 2362, ng = 7549, ng = 2395, 
nio = 8254, Tin = 1664, ni2 = 6653, nig = 3078, 
ni4 = 2739, 7115 == 5817, and 7ii6 = 1398. The acci- 
dental coincidences, which appear at the same period as 
the pump pulse period, have been subtracted. In this 
paper, all the coincidence measurement data show acci- 
dental subtracted values. 

Applying the single-ququart tomography algorithm, 
the experimentally reconstructed ququart density matrix 
is obtained to be. 



pure 
rexp 



( 0.5138 0.2749 + 0.0523i 0.3236 + 0.1308i 0.1643 -h 0.1026i \ 

0.2749 - 0.0523i 0.1590 0.1887-1- 0.0418i 0.1004 + 0.0379i 

0.3236-0.13081 0.1887- 0.0418i 0.2463 0.1259 + 0.0224i 

V 0.1643 -0.1026i 0.1004- 0.0379i 0.1259 - 0.0224i 0.0777 

I 



(6) 



We obtain Tr\p\,j,A = 0.962 which means that the ex- 
perimentally reconstructed state is not an entirely pure 
state. The entropy of the experimentally reconstructed 
density matrix is found to be S* = 0.052. 

Clearly, the experimentally prepared ququart state is 
somewhat different from what we initially intended to 
prepare and this is reflected in the state fidelity F = 
0.938. Figure [H shows the graphical representations of 



J 



the real parts of the experimental and theoretical ququart 
density matrices for this experiment. 

The errors in the experimentally reconstructed den- 
sity matrix can be estimated as follows. Assuming that 
the measured coincidences n^, follow Poissonian statis- 
tics, i.e., the ensemble average of the uncertainties satisfy 
Srii, — y/nZ^ the error in the reconstructed density matrix 
is calculated to be, 



Impure ^ 
rexv 



( 0.0066 0.0042 - 0.0034i 0.0069 + 0.0023i 0.0083 + 0.0028i \ 

0.0042 + 0.0034i 0.0036 0.0056 - 0.0035i 0.0039 + 0.0013i 

0.0069 - 0.0023i 0.0056 + 0.0035i 0.0046 0.0029 - 0.0024i 

V 0.0083 - 0.0028i 0.0039 - 0.0013i 0.0029 + 0.0024i 0.0026 / 



r 



(7) 



Mixed state ququart 



different degrees of mixedness for a biphoton ququart, we 



The double-crystal scheme shown in Fig. [T] is used to 
prepare a ququart in a completely mixed state. To study 



Re [p..,] 




Re[p„„] 




FIG. 3: Real parts of two completely mixed state ququarts 
with different entropy values. These states are generated us- 
ing the double-crystal scheme shown in Fig. [l] with differ- 
ent pump polarization, (a) For the state shown in eq. ((Ojl. 
S{pexp) = 0.394. (b) For the state shown in eq. ([T^ . 
S{pexp) = 0.504. 



studied two cases: the ququart prepared with the pump 
polarization of 30° and of 45° from the horizontal plane. 
Let us first discuss the case of 30° pump polariza- 
tion. In this case, the horizontally oriented type-I BBO is 
pumped more strongly than the vertically oriented type-I 



J 



BBO. This setting, therefore, produces an unequal mix- 
ture of the ququart basis states |3) and |0). 

The theoretical ququart density matrix for this case, 
as in eq. ([H]), takes a very simple form, 



jmix 
rtheory 



/ 0.2500 \ 





V 0.7500/ 



(8) 



The theoretical density matrix is clearly not pure since 
Tr[p^theory] = 0-625 and S = 0.406. 

For this state, the tomographic coincidence measure- 
ment resulted (for 180 s) rii = 1911, n2 = 46, na = 6287, 
m = 34, ns = 795, ng = 3562, ny = 3005, ng = 1048, 
ng ^ 1911, nio = 2061, nn = 2321, ni2 = 981, 
nu = 3141, 7114 = 3154, nig = 973,and nig = 2220. 

The experimentally reconstructed ququart density ma- 
trix is, compared with theoretical one, a bit more com- 
plex. 



/ 0.2300 

0.0024 -0.0009i 
0.0211 -I- 0.0007i 

V -0.0015 -0.0020i 



0.0024 -I- 0.0009i 

0.0057 
-0.0006 -0.0017i 
-0.0572 -h0.0019i 



0.0211 -0.0007i -0.0015 + 0.0020i\ 

-0.0006 + 0.0017i -0.0572 -0.0019i 

0.0041 0.0069 4- 0.0018i 

0.0069 -O.OOlSi 0.7571 / 



(9) 



We obtain Tr[p1 ] — 0.634 and S = 0.394 for the experimentally reconstructed ququart density matrix. The fidelity 
is calculated to be F = 0.987 and Fig. [3fa) shows the graphical representation of the real part of the experimentally 
reconstructed mixed state ququart shown in eq. ([9]). 

The experimentally reconstructed density matrix has the following inherent errors due to the fluctuations of the 
coincidence count rate. 



A mix 

rexp 



/ 0.0053 0.0027 -0.0027i 0.0030 -h 0.0023i 0.0071 + 0.0044i \ 

0.0027 + 0.0027i 0.0008 0.0044- 0.0070i 0.0044 + 0.0053i 

0.0030 -0.0023i 0.0044 + 0.0070i 0.0007 0.0048 - 0.0048i 

V 0.0071 - 0.0044i 0.0044 -0.0053i 0.0048 + 0.0048i 0.0096 / 



(10) 



r 



Let us now discuss the case in which the pump polar- 
ization is 45°. In this case, the two ququart basis states 
1 3) and |0) are equally excited. Therefore, the theoretical 
ququart density matrix is given as. 



mix 
Pth~ 



eory 



/ 0.5000 \ 





V 0.5000/ 



(11) 



The theoretical ququart density matrix results 



Tr[pfheory] = 0.5 and S = 0.5. Note that for a 
complete mixture of all four ququart basis states, 
Tr[p^] = 1/4 and 5 = 1. 

For the equal mixture of |3) and |0), the tomographic 
coincidence measurement outcomes are (for 180 s) ni = 
3442, n2 = 30, rig = 3983, n^ = 23, n^ = 1621, ng = 
2358, uy = 1950, ng = 1895, Tig = 1906, riio = 1973, 
mi = 1959, 7112 = 1840, 7113 = 2040, 7114 = 2026, 7115 = 
1809, and riig = 1909. 

The experimentally reconstructed ququart density ma- 
trix for this case is found to be, 



/ 0.4584 0.0142 + 0.0048i 0.0253 + 0.0162i 0.0158 - 0.0097i \ 

0.0142 - 0.0048i 0.0041 0.0012 + 0.0004i -0.0406 + 0.0118i 

0.0253 -0.0162i 0.0012 - 0.0004i 0.0031 0.0024 + 0.0029i 

V 0.0158 + 0.0097i -0.0406 - 0.0118i 0.0024- 0.0029i 0.5313 / 



(12) 



The experimentally reconstructed ququart density matrix is characterized by Tr[p1 ] — 0.499 and S — 0.504. The 
fidehty is calculated to be F = 0.989 and Fig. ^h) shows the graphical representation of the real part of the 
experimentally reconstructed mixed state ququart shown in eq. (fT2|) . 

The error in the experimentally reconstructed density matrix, due to the count fluctuations, is found to be. 



Ap; 



mix 
exp 



( 0.0078 0.0040 - 0.0039i 0.0043 + 0.0036i 0.0075 + 0.0047i \ 

0.0040 + 0.0039i 0.0007 0.0047- 0.0074i 0.0038 + 0.0047i 

0.0043 - 0.0036Z 0.0047 + 0.0074i 0.0006 0.0042 - 0.0042i 

V 0.0075 - 0.0047i 0.0038 - 0.0047i 0.0042 + 0.0042i 0.0084 / 



(13) 



r 



Re[p„,] 




FIG. 4: A ququart with a partially-mixed state between a 
pure state and a completely mixed state. The calculated 
fidelity is F = 0.878. We obtain Tr\pi^^ = 0.483 and 
S{p,-,-,) = 0.551. 



Partially-mixed state ququart 



A more general state of a ququart, as in eq. (|4]), can 
be prepared by transforming an initial mixed state of the 
form shown in eq. ([3]). The same double-crystal scheme 
was used for generating the SPDC photon pair and, in 
this experiment, the pump polarization was 45°. A zero- 
order half-wave plate designed at 823.5 nm was used at 
the fast-axis angle of 22.5° from the vertical axis. 



The ququart density matrix calculated for this setting 



IS, 



J 



int 
r theory 



0.2500 -.0086 0.2414 + 0.0644i \ 

0.2500 0.2414 -0.0644i 0.0086 

-0.0086 0.2414 + 0.0644i 0.2500 



\ 0.2414 - 0.0644i 



0.0086 



r 







0.2500 



(14) 



The above theoretical ququart density matrix results 
'^^[p'iheory] — ^-S and S — 0.5 which are equal to the val- 
ues for the mixed state in eq. (jlip . This means that, even 
though the ququart state has been unitarily transformed 
to a new one with the help of a wave plate, the entropy 
of the system has not been changed. To actually reduce 
the entropy of the system, it is necessary to remove the 
temporal distinguishability between the amplitudes |0) 
and 1 3) which is introduced due to ultrafast pumping of 
the BBO crystals [27|. This can be accomplished, for ex- 
ample, by inserting a properly oriented quartz plate of 



the exact thickness in the path of the pump laser. 

For the partially-mixed biphoton ququart state, the 
projection measurement resulted the following coinci- 
dence counts (for 180 s): ni = 1760, n2 = 1730, 
na = 1733, 714 = 1839, ng = 1687, ng = 1630, ny = 1758, 
rig = 1961, ng = 817, nio = 3029, nn = 1008, 
ni2 = 1701, ni3 = 1940, tim = 1944, 7115 = 1692, and 
niG = 1192. 

The experimentally reconstructed density matrix is 
found to be. 



/ 0.2493 -0.0077 -0.0007i 0.0234 + 0.0126i 0.1793 + 0.1665i \ 

i„4 _ -0.0077 + 0.0007i 0.2433 0.2032 + 0.1184i 0.0142 - 0.0036i 

Pexp- 0.0234 -0.0126i 0.2032 - 0.1184i 0.2590 0.0333- 0.0019i 

V 0.1793 -0.1665i 0.0142 + 0.0036z 0.0333 + 0.0019i 0.2542 / 

which shows Tr[pl^p] = 0.483 and S = 0.551. 

The inherent fluctuations of the count rate introduce the foUowing error in the reconstructed density matrix, 



(15) 



Art"* = 

rexp 



( 0.0059 0.0042 - 0.0042i 0.0046 + 0.0039i 0.0098 + 0.0036i \ 

0.0042 + 0.0042i 0.0059 0.0065 - 0.0054i 0.0043 + 0.0040i 

0.0046 - 0.0039i 0.0065 + 0.0054i 0.0061 0.0042 - 0.0042i 

V 0.0098 - 0.0036i 0.0043 - 0.0040i 0.0042 + 0.0042i 0.0059 

I 



(16) 



Figure [4] shows the real parts of the experimentally 
reconstructed and the theoretical density matrices. The 
calculated fidelity is somewhat low F = 0.878 in this 
case. 



DISCUSSION 



Pump 
Pulse 




Ptheory 



The fidelity F quantifies how close the experimentally 
prepared ququart state is to the one we intended to pre- 
pare, i.e., F quantifies the overlap between the theoreti- 
cal density matrix and the experimentally reconstructed 
density matrix. In section 4, we have analyzed errors 
introduced to the experimentally reconstructed density 
matrices due to the fluctuations of the count rates. As 
we have seen in eqs. (8), (11), (14), and (17), however, 
these fluctuations contribute very small errors to the re- 
constructed density matrices. 

There are a number of external experimental factors 
which could strongly affect the fidelity. First, errors in 
the angular settings of the wave plates and polarizers 
used for the projection measurement. Significant im- 
provement is possible by moving from hand-operated op- 
tic holders that are graduated in 1° ~ 2° increment to 
motorized holders. Second, less-than-ideal spatial mode 
matching for the photon pairs coming from the two dif- 
ferent crystals. Mode matching can be implemented by 
adding a spatial- filter or a short-piece of single-mode fiber 
and this should improve the fidelity substantially. Third, 
inaccurate transformation matrices of the DBS. For the 
ququart tomography discussed in the previous section, 
accurate experimental reconstruction of the ququart den- 
sity matrix requires the full knowledge of the polariza- 
tion state change induced by all optical elements. In 
our experiment, the custom-made dichroic beam split- 
ter (DBS) exhibited unexpected polarization-changing 
behaviors. It was found that the polarization states 
were changed both for the transmitted and the reflected 
beams. To account for the DBS behaviors, we carried out 



FIG. 5: Scheme for preparing a general biphoton ququart. 
Pump HWP is used to control the relative magnitude of the 
SPDC amplitudes generated at the BBO. Four BBO (two 
type-I and two type-II) crystals are needed for this scheme. 



the Stokes parameter measurements for the transmitted 
and the reflected beams for six different input polariza- 
tion states (|7J),|y),|45°),|135°),|i?), and \L)). From 
these measurements, it was possible to deduce the 2x2 
DBS transformation matrices for the transmitted and the 
reflected modes. The experimental ququart density ma- 
trices shown in the previous section were reconstructed 
using the DBS matrices and therefore experimental er- 
rors introduced to the DBS matrices should have slightly 
affected the fidelity. 

In our experiment, preparation of general states of 
a ququart was demonstrated using the double-crystal 
scheme. This scheme, although easy to setup, has poten- 
tial difficulties in complete control of the ququart states. 
For example, independent decoherence control for ampli- 
tudes |1) {\Hx,,Vx^)) and |2) {\Vx,,Hx^)) using a bire- 
fringent medium is difficult as both amplitudes contain 
horizontal and vertical polarization components. 

For a complete control of the ququart state, i.e., to 
prepare a ququart state with arbitrary values of Trlp^] 
and S, we can envision a four-crystal scheme in which a 
SPDC source is placed in each arm of a four-path Mach- 
Zehnder interferometer as shown in Fig. [51 This scheme 
requires four crystals (two type-I SPDC sources and two 
type-II SPDC sources) and the four-path interferometer 



must be made stable for an accurate phase control. De- 
coherence control can be accomplished by controlling the 
effective beam paths of the interferometer arms using the 
tunable delays installed at each of the beam path. (Note 
that arbitrary pure ququart states can be prepared using 
just two crystals, see Ref. 20|,|2l|.) 



It is, however, possible to design experimental schemes, 
for preparing arbitrary (mixed, partially-mixed, and 
pure) ququarts states, which are less complicated than 
the four-path intcrferometric scheme shown in Fig. [5] In 



the following, we discuss two such experimental schemes, 
each of which are suited for particular ququart states in 
need. 



Scheme based on Mach-Zehnder interferometer 

The basic idea for the new ququart preparation scheme 
is based on the observation that the pure ququart state 
in eq. ([2|) can be re- written as. 



I 



(17) 



where (pij {i,j — 0, 1, 2, 3) is the relative phase for the i 
and j ququart basis states. Here, we have put together 
the two ququart amplitudes that can be prepared with 
type-I SPDC (the first two terms) and the other two 
ququart amplitudes that can be prepared with type-II 
SPDC (the last two terms). 

Experimentally, the first two terms in eq. (|17p can 
be prepared with two orthogonally oriented type-I BBO 
crystals placed in tandem as in Fig. [1] The second two 
terms can then be prepared with two similarly placed 
type-II BBO crystals. Therefore, coherently or incoher- 
ently combining these two experimental schemes will al- 
low us to prepare an arbitrary ququart states. 

The experimental scheme realizing this idea is shown 
in Fig. [6] A Glan-Tompson prism (GP), transmitting 
the horizontally polarized component of the UV pump 
and reflecting the vertically polarized component, serves 
as the input beam splitter for the Mach-Zehnder inter- 
ferometer. The reflected pump laser, after passing the 
birefringent compensator (Comp.) and a half-wave plate 
(HWP2), pumps a set of two orthogonally oriented type- 
I BBO crystals. A quartz plate (QPl) compensates the 
group velocity delay between the ordinary and the extra- 
ordinary polarized photons emitted from the pair of BBO 
crystals. The residual pump laser is removed by a uv mir- 
ror (UVM) and two quartz plates QP2 can be tilted along 
their optical axes to introduce a phase shift (^os) be- 
tween horizontally and vertically polarized type-I bipho- 
tons 27| . The dichroic mirror (DM) is designed to trans- 



HWPl GP Comp. HWP3 QP3 

+- B ■ I — V/-\ 



mit the biphoton wavelengths but to reflect the uv pump 
beam which comes from the upper path of the Mach- 
Zehnder interferometer. 

The uv pump beam in the upper path of the Mach- 
Zehnder interferometer goes through a birefringent com- 
pensator, phase plates (QP3) for adjusting the phase 
((/>i2) between horizontal and vertical components of the 
pump beam, and a piezoelectric translator (PZT) which 



Comp. 




Type-I 
BBO QPl UVM QP2 



Type-II 
BBO QP4 UVM 



I DDD I \/ ^K^^+ 



Pthe 



FIG. 6: Proposed setup for arbitrary ququart preparation. 
This scheme reduces the four-path interferometer shown in 
Fig. E] to a Mach-Zehnder interferometer, making the exper- 
imental implementation much easier. QP2, QP3, and PZT 
adjust the phase terms in eq. (|17p . (^03, 012, and (^oi, respec- 
tively. 



introduces a relative phase shift {(t>oi) between the up- 
per and the lower paths of the interferometer. The uv 
pump from the upper path, upon reflection at the DM, 
serves as the pump for the set of two orthogonally ori- 
ented type-II BBO crystals. The group velocity delay 
between biphoton amplitudes from the first and the sec- 
ond type-II BBO crystals are then compensated by the 
crystal compensator (QP4). 

At the output of the experimental setup shown in 
Fig. [SI an arbitrary ququart states based on the biphoton 
polarization states of frequency-nondegenerate SPDC is 
prepared. We note that the scheme is loosely based 
on the biphoton qutrit setup demonstrated in Ref. [18|, 
where three SPDC crystals are exploited for generation 
of frequency-degenerate coUinear biphotons. 



Scheme based on frequency non-degenerate and 
non-collinear regime of SPDC 

The interferometric scheme proposed in the previous 
section, although straightforward, might not be practi- 
cal as the scheme inherits high phase sensitivity of the 
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Type-I 
WP QP BBO 




FIG. 7: Proposed experimental setup for preparation of an 
arbitrary pure ququart state. DBS represents the dichroic 
beam splitter. 



Mach-Zehnder interferometer which is not desirable for 
the purpose of preparing ququarts. Is it then possible 
to prepare an arbitrary ququart states using less compli- 
cated experimental scheme? The answer to this question 
is found to be positive, at least for the pure ququart states 
by using the biphoton polarization entangled states with 
controllable two-qubit concurrence C . 

Let us first discuss the problem in a rather abstract 
form. In some arbitrary chosen four-dimensional com- 
putational basis, an arbitrary pure ququart state can be 
written as eq. ((2]). Since the biphoton ququart is in fact 
formed with two polarization qubits as defined in eq. ^ , 
we can make use of the fact that a unique set of or- 
thonormal states of the two subsystems (two polariza- 
tion qubits) \Ai) and \Bi) (i — 1,2) exists such that the 
biphoton pure state in eq. ^ can be expressed in the 
form 



m ^ VxI\Ai)\A2) + ^^\B,)\B2). 



(18) 



This is known as the Schmidt decomposition. The 
Schmidt coefficients Xi ^re eigenvalues of the reduced 
density matrices of the subsystems. If a biphoton 
ququart state can be decomposed in this way (although 
the Schmidt coefficients and the Schmidt basis states may 
be different), one should be able to prepare an arbitrary 
pure ququart state, provided that the Schmidt coeffi- 
cients and the Schmidt basis are experimentally control- 
lable. The control of the Schmidt coefficients requires 
non-local unitary transformations which affect both sub- 
systems but it is possible to use only local operations to 
switch between the Schmidt basis. For biphoton ququarts 
we discuss in this paper, these operations turn out to be 
rather simple as the subsystems are polarization qubits. 
The proposed setup to implement this idea is shown 
in FiglTl A set of two orthogonally oriented type-I BBO 
crystals, cut for frequency-nondegenerate non-collinear 
type-I SPDC, are pumped by a uv laser whose polar- 
ization is controlled by a half- wave plate (WP). A set 
of quartz plates (QP) introduces the relative phase be- 
tween the horizontal and the vertical components of the 
uv pump. The initial polarization state of the biphoton 
generated in this process can be written as 



where the relative magnitude of the coefficients is con- 
trolled by WP and the relative phase is controlled by 
tilting QP. For this biphoton state, the state of individ- 
ual subsystem (qubit) is then given as. 



P,^ximm+X2m{v,\, 



(20) 



where the subsystem index i — 1,2. 

So far we have produced a state with a given set of 
Schmidt coefficients (i.e., a given concurrence), but in 
a fixed (H-V) basis. The use of a quarter-wave plate 
(QWP) and and a half-wave plate (HWP) on each of the 
photons then locally and unitarily transforms the polar- 
ization states of the individual photon. Since the uni- 
tary transformation conserves the inner product, there 
always exists a state orthogonal to a given one in the two- 
dimensional Hilbert space. In other words, in eq. (|20p . if 

\H,) ^ |A,), then \Vi) ^ \B,) with {A,\B,) - 0. As a 
result, the state eq. (fT9)) becomes 



(7i®C/2 



xI\Ai)\A2) + ^X^\Bi)\B2 



(21) 



where C/i, for example, refers to the unitary polarization 
transformation for photon 1. Equation (|2ip represents a 
general form of a pure biphoton polarization state or a 
pure biphoton ququart state. 

The scheme shown in Fig.[7]is interferometrically more 
stable than the one in Fig. [5] due to the fact that the 
photon pair from each type-I BBO crystal goes through 
the same optical paths. In principle, to prepare the state 
in eq. (PT|) , the two paths in Fig. [7] should only be equal 
up to the coherence length of the pump laser. If the pump 
is broadband, however, it becomes necessary to further 
erase the temporal distinguishability between amplitudes 
from the first and the second BBO crystals, for example, 
by using a set of compensating crystals 27 1. 



CONCLUSION 

We have demonstrated that an arbitrary general single- 
ququart state can be prepared in a simple and con- 
trollable way by using the biphoton polarization state 
of ultrafast-pumped coUinear frequency-nondegenerate 
SPDC. In addition, we have proposed two additional 
schemes which can be applied for arbitrary (pure and 
mixed) ququart state preparation. 

Compared to other multi-dimensional quantum sys- 
tems, the biphoton ququart is easier to prepare and char- 
acterize and states other than pure states can be prepared 
easily. Furthermore, it is possible to prepare a multi- 

We, 



22, m 



Ti\Hx,)\Hy, 



X^\Vx,)\Vx,), (19) 



ququart entangled state linear optically 
therefore, believe that the general ququart state prepara- 
tion scheme analyzed in this paper will find applications 
in quantum key distribution and quantum information 
processing. 
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